Voltage and temperature dependencies of conductivity in gated graphene 
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The resistivity of gated graphene is studied taking into account electron and hole scattering by 
short- and long-range structural imperfections (the characteristics of disorder were taken from the 
scanning tunneling microscopy data) and by acoustic phonons. The calculations are based on the 
quasiclassical kinetic equation with the normalization condition fixed by surface charge. The gate- 
voltage and temperature effects on the resistance peak, which is centered at the point of intrinsic 
conductivity, are found to be in agreement with the transport measurements. 
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Numerous features of electronic properties of graphene 
(see discussion and Refs. in 1 and 2) are caused by both a 
neutrino- like dynamics of carriers [3( , which is described 
by the Weyl- Wallace model [3], and a substantial mod- 
ification of scattering processes. For instance, the gap- 
less energy spectrum and specific character of scattering 
are responsible for the peak of resistance of the gated 
graphene sheet (Fig. la). Such a peak was observed in 
Refs. 1 and 5-9 and one appears due to the transforma- 
tion between n- and p-types of conductivity through the 
intrinsic region with the chemical potential in the vicinity 
of the band cross-point. Moreover anomalous tempera- 
ture dependence was found: the peak resistance decreases 
with temperature while the tail resistance increases. To 
the best of our knowledge, a qualitative description of 
the gate-voltage and temperature effects on conductiv- 
ity is not performed yet. Much attention was attached 
to the minimal conductivity phenomena, see last papers, 
[To| and the metallic regime of conductivity was also ex- 
amined. pH Recently, the near-maximum shape of the 
graphene resistance was considered assuming the scat- 
tering by a remote charge impurities layer is a dominant 
scattering mechanism. Q Since the high-density layer is 
not likely to be present, an additional scattering due to 
disorder was suggested in Ref. 12, but the only metallic 
regime of conductivity was discussed. 

In this Brief Report, we calculate the gate- voltage and 
temperature dependencies of the resistivity peak tak- 
ing into account structural inhomogeneities of graphchc 
[both long-range inhomogeneities and point defects, Figs. 
lb and lc, respectively] and the acoustic phonon scat- 
tering. Our model is justified by the recent scanning 
tunneling microscopy (STM) measurements, [l3| where 
both the imperfections with the lateral scale of 5—9 nm 
and point defects were reported. The long-range disorder 
(DL) is described by the potential t/ x with the Gaussian 

correlation function (U X U X '} = Ui exp{ — [(x — x.')/l c ] 2 } 
where Ui is the averaged energy and l c is the correlation 



length. The short-range defects (DS) of sheet density rid 
are approximated by the potential t/dA(x) where A(x) 
is a function localized over the scale I 2 ,. The main contri- 
bution to the acoustic phonon scattering appears due to 
the deformation interaction with longitudinal vibrations, 
[Til [l5| DV • u x , where D is the deformation potential 
and u x is the displacement vector of LA-mode. 
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FIG. 1: (a) Schematic plot of the resistance at zero (solid) and 
finite (dotted) temperatures versus gate voltage. Nonideal 
sheet of graphene (6) with a long-range inhomogeneities (gray) 
and (c) with point defects marked by balls. 



The peak of resistance is originated from the distinc- 
tions of scattering caused by DS (or LA) and DL mech- 
anisms. Since a short-range relaxation rate (due to DS 
contribution or LA phonon scattering at high tempera- 
tures) is proportional to the density of states [ll[, the 
resistance R does not change with concentration or gate 
voltage V g . Besides, the LA contribution is proportional 
to temperature of phonons T, and R increases with T. 
In contrast, the DL scattering becomes suppressed if an 
electron (hole) momentum exceeds h/l c . Due to this, 
R decreases when V g or T (i.e. energy of carriers) in- 
creases. Thus, a combination of DL- and DS- (LA-) 
mechanisms allows an explanation of an experimental be- 
haviour shown in Fig. la. 

Scattering between electron or hole states [l6[ with 
two-dimensional momenta p and p' is described by the 
transition probabilities W^} where j = DL, DS, and 
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LA. Within the Born approximation, W^^, is wril Ion in 



2 



the golden rule form, 
2tt 



K3 = -ir w l U) *^ w {p-p') 



(i) 



where the overlap factor = (1+cos 9)/2 with 9 = p', p 
describes the suppression of backscattering processes, 
[13] vqw — 10 s cm/s is the characteristic velocity of the 
linear dispersion laws, iv^wP, with the upper and lower 
signs corresponding to electrons and holes, respectively. 
The matrix element Wq , which depends on the mo- 
mentum transfer ftq = p — p', is given by 



; -(«/c/2r 

2T/hu> q 



l<V a 



(DL) 
(DS) 
(LA). 



(2) 



Here the DL and DS mechanisms of scattering by static 
disorder are described by the Fourier transformations of 
{U X U X ') and [7dA(x). The LA probability is expressed 
via the matrix element \C q \ 2 = D 2 TiLo q / (2p s s 2 L 2 ), where 
s is the sound velocity, uj q = sq is the phonon frequency, 
p s is the sheet density of graphene, and L 2 is the nor- 
malization area. Within the equipartition condition, at 
foujq <C T, we replace the Planck distribution of phonons 
by T/fiLO q . 

Under the steady-state electric field, E, the electron 
(e) and hole (ft) distribution functions, f ep and fh P , are 
governed by the quasiclassical kinetic equations (here and 
below k = e, ft), 



dfkp 
dp 



(3) 



with the standard collision integral for the elastic scatter- 
ing. Using the electron and hole velocities ±vswP/p, 
one obtains the following expression for the current den- 
sity: 



I 



Aevgw \ - P 
L 2 ^ p 



{fep ~l~ fhp^) 



(4) 



where factor 4 appears due to the spin and valley degen- 
eracy. The distribution functions f ep and fh p are also 
related to each other and to the sheet charge density, 

Qs, 



4e 



(5) 



Note, that Q s can be expressed via the gate voltage as 
V g = eQ s /ATTd. Below we assume the SiC>2 substrate of 
the width d =3 nm with the dielectric permittivity e ~3. 

Within the linear approximation, the solution of Eq. 
(3) can be seached in the form fu p = F^p + A/fc p , where 
Fk P = {exp[(wi3wp T n)/T} + l} -1 are the equilibrium 
Fermi distributions with the chemical potential p. The 
asymmetric parts of distribution functions, A/^p, can be 
obtained from Eq. (3) in the form 



AA P = - 



(eE-p) 

P 



JB(m) 



dFi 



kp 



dp 



(6) 



with the momentum relaxation time, Tp^ m \ is given by 

-l 



the sum of relaxation rates r, 



fi(m) 



E 7 (i/ 



where 1/t$ U) = £\ , W p ^>,(l - cos 9). Substituting Eq. 



Mi) 



p,p 



(6) into Eq. (4) one obtains the resistivity, which is in- 
troduced by I = E/i?, as follows: 



1 

R 



e 2 vew 
irh 2 



dp v** £ 



dFi 



hp 



dp 



(7) 



and the momentum relaxation time here do not depend 
on k. It is convenient to present the relaxation rates as 



1 



„ B (J) 



p / * (plc/h) (DL) 

1 (DS, LA), 



where the averaging over angle gives the function 



z 2 /2 

~ h \ 2 



(8) 



(9) 



with the first-order Bessel function of an imaginary ar- 
gument, h(x). The characteristic velocities, v^J\ in Eq. 
(8) are introduced as follows: 



Mi) 



so that Vsi^ <C vqw for the parameters used below. 

Next, we simplify the expression (7) using the di- 
mensionless resistivity, Re 2 /irfi. We consider the case 
of degenerate carriers, \p\ 3> T, when Eq. (5) gives 
the square-root dependence of p on the gate-voltage: 
p = hvftw V 'eV g I A\e\d. Since J2k(~dFkp/dp) can be re- 
placed by 5- function, Eqs. (7)- (9) yield the dimensionless 
resistivity in the following form: 



1 



4ft v B w 



(DL) 
(DS) 

D 2 T/( Ps s 2 ) (LA), 



(10) 



R 



(11) 



with the temperature dependent LA contribution 
v% L /vsw- Another simple expression can be obtained 
for the case of the intrinsic regime of transport, Q s = 0, 
when p—0 and Fk p — + [cxp(vuwp/T) + 1] _1 in Eq.(7). As 
a result, the inverse sheet resistance is expressed through 
the dimensionless integral, 



nil 
lie 2 

TL 



v 



S(DL) 



dx 



v SiW h 



VBW JO 

„M DS ) 

Use. 



cosh a; 
-i -l 



(12) 



VQW 



Here, not only the ratio Vsi L /«ew increases with the 
temperature but also the DL contribution becomes sup- 
pressed, if Tl c /vaw ft > 1- 



S(LA 



VQW 



FIG. 2: (Color online) Normalized resistance R/R m versus 
dimensionless gate voltage V g /V s plotted for parameters of 
ref. 1 (dotted curve, hexagons), Ref. 7 (solid curve, triangles) 
and Ref. 8 (dashed curve, squares). 



We turn now to discussion of the experimental data 
using the parameters of Eq. (2) in order to fit the gate- 
voltage and temperature dependencies of resistance ob- 
tained in Refs. 1, 7, and 8 (similar low-temperature re- 
sults can be obtained from Refs. 6 and 9). Under the 
description of the LA scattering we use the known pa- 
rameters s ~ 7.3 • 10 5 cm/s and p s ~ 7 • 10~ 8 g/cm 2 . 
In order to obtain the temperature-dependent contribu- 
tion to the tails of peak [about 100 in Fig. 2a of Ref. 

7 and in Fig.l of Ref.8] one needs the deformation po- 
tential D ~12 eV which is about 75% of the graphite 
value. [13, [TtJ For such a set of parameters, one obtains 
Vs { c LA) ~ 8 • 10 5 cm/s at T=300 K from Eq. (11). Since 
ty(z — > oo ) = 0, the temperature independent part of 
the tails, which is about 400 fl in Fig. 1 of Refs. 1 and 

8 and about 200 fi in Fig. 2a of Ref. 7, appears due 
to the DS contribution with the characteristic velocities 
Vs { c DS) ~ 3.2 • 10 6 cm/s and 1.6 • 10 6 cm/s for Refs. 1 and 
8 and Ref. 7, respectively. 

In order to fit the gate-bias dependencies, in Fig. 2 
we plot the normalized resistance R/R rn versus the di- 
mensionless gate-bias, V g /V s , where R rn = Trh(Vsi DL ' + 
Vsi DS * > )/vftwe- 2 and V s — 4|e|d/e£j?. Here we used the 
values of Vsc obtained above and the 5 nm correla- 
tion length, which is taken from the STM measurements 
[l3| , for all the cases. [3, 0, A good agreement is ob- 
tained for Vsi DL ^ ~ 3.2 • 10 7 cm/s (Refs. 7 and 8) and 
5 ■ 10 7 cm/s. Q Thus, we have found the scattering pa- 
rameters of the relaxation rates given by Eqs. (8) and (9) 
for three different samples and one can plot Vp versus 
energy v$wP (see Fig. 3). The velocities Vp^ DL DS ^ are 
connected to the characteristic potentials Ui and Ud as 
follows. For the long-range disorder, we obtain U\ ^100 
meV or 80 meV for Ref. 1 or Refs. 7 and 8, respectively, 
i.e. U i is less than the interlayer coupling energy ^350 
meV. The energy Ud can be estimated as ~3 eV (of the 
order of the intralayer coupling energy), if we use Id —2 
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FIG. 3: Relaxation rates l/r^ versus energy, vswP- (1) 
j — DL for parameters of Refs. 7 and 8, (2) j — DL for 
Ref.l, (3) j = DS for Refs. 1 and 7, (4) j = DS for Ref. 8, 
and (5) j = LA for Refs. 1, 7, and 8. 



nm and ndl\ ~ !%■ The suppression of resistance with 
V g takes place if the chemical potential exceeds the 100 
meV threshold, see Fig. 3 (p. ^50 meV at V^=10 V for 
the structures under consideration). Thus, both voltage- 
induced and temperature-induced (below 700—1000 K 
where the optical phonon contribution becomes essen- 
tial) suppression of resistance takes place. 
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FIG. 4: (Color online) Temperature dependence of resistance 
R/Rm for the intrinsic regime of conductivity, fj, = 0. The 
experimental data taken from Refs. 7 and 8 are marked as in 
Fig. 2. 

Due to lack of data on the temperature dependency 
of resistance, we only discuss here the suppression of the 
resistance peak value when temperature increases. We 
use Eq. (12) which governs the intrinsic regime of con- 
ductivity and the dependency R(T) is in agreement with 
the room temperature data of Refs. 7 and 8 without 
any additional parameters. Figure 4 shows an essential 
suppression of R when the ratio Tl c /vsw^- exceeds unit 
(at T >300 K for the parameters used). An extrapola- 
tion for the high-temperature region shows an essential 
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suppression of the DL scattering at T >1000 K. Thus, 
a minimum of resistance, which is caused by such a sup- 
pression and an optical phonon contribution, should take 
place for the temperature region above ~ 500 K. 

Let us discuss the approximations used in our cal- 
culations. The main assumption of our model is the 
phenomenological description of the elastic scattering by 
long- and short-range structural inhomogeneities. While 
the in-plane scales of disorder, l c and Id, are clear and jus- 
tified by Ref. 13, a nature of inhomogeneities and a value 
of the phenomenological energies [Ui and Ud in Eq.(2)] re- 
quire an additional analysis. Other approximations made 
are rather standard and generally accepted. We have ne- 
glected the screening effects because it is not a leading 
effect for a spatial ranges up to 5 nm. The only scatter- 
ing by LA mode is taken into account because another 
contributions are weaker (see Refs. 14 and 17 and refer- 
ences therein) . Since Vsi *C Vbw, the (anti)localization 
phenomena can be neglected if Re 2 /nh <C 1 (but it 
may be essential in Ref. 5 and in Fig. 3 of of Ref. 7, 
where R m ~ 10 Kfi). In addition, the quantum dynam- 
ics restrictions, like a Zitterbewegung effect [20j or the 
field-induced interband tunneling, are only essential in 
the vicinity of the cross-point. Since R is weakly changed 
near this point, one can use the results obtained as phe- 
nomenological dependencies. However, a justification for 
the scales ~1 V or tens of kelvins is beyond of the con- 
sideration performed. 

In summary, we have found that the main factor, which 
forms the resistance peak centered at the intrinsic con- 



ductivity region, is a structure disorder (e.g., defect clus- 
ters, bilayer islands, or a substrate disorder [5lJ) with 
the scale around 5 nm for the samples analyzed. Short- 
range defect contribution and acoustic phonon scatter- 
ing are only essential for the heavily doped case on 
the tails of peak. Both gate-voltage and temperature- 
induced quenchings of the resistance peak are in agree- 
ment with the transport measurements [l], 0, Q and 
with the STM mapping of nonideal graphene sheet. [l3| 
A complete verification of the mechanism suggested re- 
quires both microscopic calculations (in order to estimate 
the fitting parameters Ui and Ud) and a further investi- 
gation of the transport phenomena. Particularly, high- 
temperature conductivity measurements, which can con- 
firm a suppression of resistance peak, as well as treat- 
ment of magnetotransport and high-frequency responses 
are necessary. 

To conclude, the results presented permit to analyze 
linear characteristics of the graphene-based field-effect 
transistor. Recent studies of such a device, see Ref. 22 
and references therein, demonstrated an essential gate 
modulation. However, the data available are not suffi- 
cient to determine the scattering mechanisms. We be- 
lieve that t result obtained will stimulate further investi- 
gations in order to understand transport phenomena and 
to improve device characteristics. 
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the Japan Science and Technology Agency, CREST. 
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